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Using a random-matrix approach and Monte-Carlo simulations, we generate scattering matrices 
and cross sections for compound-nucleus reactions. In the absence of direct reactions we compare the 
average cross sections with the analytic solution given by the Gaussian Orthogonal Ensemble (GOE) 
triple integral, and with predictions of statistical approaches such as the ones due to Moldauer, 
to Hofmann, Richert, Tepel, and Weidenmuller, and to Kawai, Kerman, and McVoy. We find 
perfect agreement with the GOE triple integral and display the limits of validity of the latter 
approaches. We establish a criterion for the width of the energy-averaging interval such that the 
relative difference between the ensemble-averaged and the energy-averaged scattering matrices lies 
below a given bound. Direct reactions are simulated in terms of an energy-independent background 
matrix. In that case, cross sections averaged over the ensemble of Monte-Garlo simulations fully agree 
with results from the Engelbrecht-Weidenmuller transformation. The limits of other approximate 
approaches are displayed. 
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I. INTRODUCTION 

For medium-weight and heavy target nuclei, nuclear 
reactions represent a very complex phenomenon because 
the number of degrees of freedom grows rapidly with mass 
number A. That fact has naturally led to the develop¬ 
ment of a statistical approach. Central to the approach 
are the concept of a fully equilibrated compound nucleus 
and the Bohr hypothesis [T] , which states that a particle 
incident on a medium-weight or heavy nucleus shares its 
energy with the target nucleons. The resulting compound 
nucleus attains statistical equilibrium, and the modes of 
decay of the equilibrated system are independent of the 
mode of formation. The postulated independence im¬ 
plies a factorization of the energy-averaged compound- 
nucleus cross section [2]. The factorization hypothesis 
holds very well at sufficiently large bombarding energies 
(i.e., in the Ericson regime) but not for isolated or weakly 
overlapping compound-nucleus resonances [314]. In that 
regime the average cross section given by the factoriza¬ 
tion hypothesis must be corrected by a “width fluctuation 
correction” factor (WFC). The WFC factor basically ac¬ 
counts for an enhancement of the elastic average cross 
section. 

In the 1970s numerous efforts were undertaken to de¬ 
rive the WFC factor [iHID] or to generate a suitable 
parametrization of the WFC factor with the help of the 
Monte-Carlo (MC) technique [S] ITIHTB] . All of these 
were guided by random-matrix theory (RMT). Inspired 
by Bohr’s idea, Wigner had introduced RMT into nu¬ 
clear physics as a means to cope with the complexities of 
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the compound nucleus (see Brody et al. [H]). In RMT, 
the nuclear Hamiltonian is assumed to be a member of 
the Gaussian Orthogonal Ensemble (GOE) of random 
matrices. Wigner himself never went as far as formulat¬ 
ing a statistical theory of nuclear reactions in terms of 
the GOE. Lacking such a theory, the above-mentioned 
approaches used approximations that were not fully con¬ 
trolled. Only in 1985 an exact closed-form expression 
for the average S matrix and for the S matrix corre¬ 
lation function based upon a GOE scattering approach 
was derived [T^, based upon the shell-model approach 
to nuclear reactions m and valid in the limit of a large 
number of resonances. In that work the S matrix is writ¬ 
ten in terms of the GOE Hamiltonian Averages 

are performed directly over the Gaussian-distributed el¬ 
ements of 

The exact results of the GOE scattering approach [15] 
apply for all values of the parameters (number of open 
channels, isolated or overlapping resonances) character¬ 
izing compound-nucleus reactions. More generally, that 
work describes universal features of quantum-chaotic 
scattering [HI ng and is, therefore, relevant also be¬ 
yond the confines of nuclear physics. However, the ex¬ 
act expression for the S'-matrix correlation function m 
involves a triple integral. The computational cost of 
evaluating that integral is quite heavy especially when 
many channels are open. That is why only few numer¬ 
ical studies have been performed in the past. Erdhner 
m and Igarasi [20] independently compared the GOE 
triple integral results with Moldauer’s method, and ob¬ 
tained good agreement. Hilaire, Lagrange, and Koning 
[2T] extended the numerical study, and applied it to some 
realistic cases where neutron radiative capture and fis¬ 
sion channels are involved. Updated parametrizations of 
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Moldauer’s method based on the GOE triple integral cal¬ 
culation are also available [22l [23] , which are of practical 
use for cross-section calculations. 

In the present paper we present a thorough analysis 
of the results of the GOE approach and a comparison 
with other, approximate methods, with the aim to un¬ 
derstand their applicability and limitation. The work is 
based upon a Monte-Garlo approach. We generate an 
ensemble of scattering matrices or cross sections. This 
is done by drawing at random the elements of 
and using these to generate the elements of the scatter¬ 
ing matrix. In this way we are able to avoid some phe¬ 
nomenological assumptions made in the past concerning 
the distribution of the decay amplitudes or of levels. In 
that respect our MG approach also differs from the one 
used by Moldauer or Hofmann et al. We average over 
the ensemble of realizations generated by the MG method 
and compare these with predictions of the exact GOE ap¬ 
proach and of other approximate methods. We are able 
to answer some important and long-standing questions 
concerning compound nuclear reactions, such as the dif¬ 
ference between energy and ensemble averages, the role 
of direct channels, the existence of correlations between 
distributions of levels and decay amplitudes, and the be¬ 
havior of the cross section in the limit of weak absorption. 


II. THEORY OF STOCHASTIC SCATTERING 


A. Compound-nucleus cross section 


The cross section for a reaction from channel a to chan¬ 
nel b is written as 


^ab — 1^2 | 


( 1 ) 


Here ka is the wave number for channel a, Qa is the spin 
factor, and the element Sab of the scattering matrix S 
consists of an energy-averaged part (Sab) and a fluctu¬ 
ating part S'^j. The energy-averaged cross section also 
consists of two parts, 


(0-o&) = -^9a{\Sab - Sabl"^) 

= ^9a{\Sab-{Sab)\^ + {\S^b\^)} ■ ( 2 ) 


The term containing |5at, — (S'a&)P describes shape elas¬ 
tic (a = b) or shape inelastic (a ^ b) scattering. The 
term containing is the average compound-nucleus 

(GN) cross section 


= ^9a{\Sab\^) . (3) 

In the first part of the paper we confine ourselves to cases 
where the average S matrix is diagonal, (Sat) = bab{Saa)- 
Then {Saa) and the shape-elastic cross section 
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SE 

aa 


^9a\l-{Saa)\^ 


(4) 


are given by the optical model. It is the aim of various 
theories of GN reactions to express the CN cross section 
in terms of {Saa) and of the transmission coefficients 

T; = l-|(^aa)|" , 0<r„<l. (5) 


These measure the unitary deficit of (S) and, thus, the 
probability of CN formation. In the second part of the 
paper we address the case when (Sab) is not diagonal. 

Bohr’s idea of the independence of formation and decay 
of the CN led to the Hauser-Feshbach formula [5] for the 
CN cross section. 


HF _ TaTb 


( 6 ) 


Corrections to that formula are conveniently expressed 
in terms of the “width fluctuation correction” (WFC) 
factor [21] . 


a 


CN 

ab 


_ TT TaTb 

~ kf^EaTc 


Wab ■ 


(7) 


Rigorously speaking, Wab should be separated into two 
parts, the “elastic enhancement factor” and the proper 
“width fluctuation correction factor” [Ij. However, for 
the comparison of various approaches it is more conve¬ 
nient to adopt the suggestion of Hilaire, Lagrange, and 
Koning and to define the width fluctuation factor 
Wab as the ratio cr^b /^ab- 

In what follows we compare several approaches to the 
calculation of and/or of the WFC factor. In chrono¬ 
logical order, these are the approach of Kawai, Kerman, 
and MacVoy [6| (KKM), the parametrization by Hof¬ 
mann, Richert, Tepel, and Weidenmiiller |5] (HRTW), 
Moldauer’s parametrization [Tl], the GOE approach by 
Verbaarschot, Weidenmiiller, and Zirnbauer m, the 
parametrization by Ernebjerg and Herman [22| . and that 
by Kawano and Talou [23] . These are briefly summarized 
in the Appendix. Hereafter we drop the kinetic and spin 
factors TTga/ka, so that all the cross sections are dimen¬ 
sionless. 

All these approaches use GOE-inspired statistical as¬ 
sumptions on the parameters of the CN resonances. In 
our comparison we use the results of Ref. m as a bench¬ 
mark. We do so because the work of Ref. [15] is the only 
one that, starting from a random-matrix model for the 
Hamiltonian of the CN resonances and using controlled 
approximations, obtains an analytical expression for 
that is valid in all regimes — from the regime of isolated 
resonances to that of strongly overlapping resonances. 


B. S matrix, K matrix, R matrix 

In order to display the connection between various 
theories of resonance reactions we recall here briefly 
the derivation of a universal expression for the S ma¬ 
trix PITB]. Specialization of that expression then yields 
the formulas used in various approaches. 
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Given a time-reversal-invariant Hamiltonian H we use 
Feshbach’s projection operators P and Q = 1 — P (where 
P projects onto all open channels labelled a,b,...) to 
write the Schrodinger equation {E — H)'i> = 0 for the 
scattering wave function 'I' in the form of the coupled 
equations 


{E - Hpp)P'i' = HpqQ'S , (8) 

{E - Hqq)Q^> = HqpP^ . (9) 

We use the standard notation, iJpp = PHP, Hpq = 
PHQ, etc. With the P space scattering wave function 
defined by 


{E - Hpp)i;i+^ = 0 , (10) 


the unitary and symmetric S matrix is given by 


Sat = - 2m 

Here is a unitary background scattering matrix de¬ 
fined by the asymptotic form of the solutions and 

P-QQ is the effective Hamiltonian in Q space, 


Pqq = Hqq + Hqp 


E+ - Hpp 


~H 


PQ 


( 12 ) 


To be useful Eqs. (11) and (12) must be specialized fur¬ 
ther. The ^-matrix approach of Ref. m and the expres¬ 
sions for S in terms of the K matrix and the R matrix use 
different such specializations. Common to these is the as¬ 
sumption that the unitary background scattering matrix 
^(0) 

is diagonal, = (5a& exp{2f())a}. We assume that 
the phases /pa are removed by the transformation Sat —t 
exp{-f(/)a}S'at,exp{-z^f,}. Then exp{-f(/)a} = V’a is 
real, and Eq. 0 becomes 

Sab = Sab - 27rf (^^o.\HpQ ^ . (13) 


For the S'-matrix approach of Ref. [15] we introduce 
an arbitrary orthonormal basis of states labeled /i in Q 
space and write 


W^a = {f^\HQp\i^a) = Wa^ = , (14) 



1 

E+ - Hpp 


Hpq\i 


(15) 


assumption throughout. With these definitions, Eqs. 0 
and (12) take the form 


Sab = Sab - 2i7rJ2Wa>.iD-^)^.W,b (17) 


where 


D - ESfj^iy 


■ I'K 




(18) 


For the iF-matrix parametrization of S we use Eq. (12) 
to define the eigenvalues E^ and eigenvectors of the 
bound compound system. 


{Ea - HQQ)Xa = 0 . 


(19) 


The states X^ produce the CN resonances in the scatter¬ 
ing process. These states correspond to a special choice 
of the basis of states /r used in Eqs. (141-( 16). The partial 
decay amplitude of state X^ into channel a is 

laa=lla = {Xa\HQp\i^a) . ( 20 ) 

Under neglect of the shift matrix A the S matrix of 
Eq. (13) can be written as 


Sab = 


l-iK 
1 -I- iK 


( 21 ) 


ab 


where 




Kab{E) = nYiY 


( 22 ) 


The R matrix is obtained by a non-standard choice of 
the projection operators P and Q. In every channel c 
(open and closed) a radius Vc is defined. The set {rd of 
all channel radii separates the internal and the external 
regions of configuration space. The operator Q projects 
onto the internal region, and P = 1 — Q. At the chan¬ 
nel surfaces of the internal region, self-adjoint boundary 
conditions are introduced with real boundary condition 
parameters Be- These define a Hermitian Hamiltonian 
Hqq and associated internal eigenvalues and or- 

thonormal eigenfunctions X^ ’ ■ The reduced width am- 
plitude ^ac is the projection of the eigenfunction X^ ' 
onto the surface of channel c, and the R matrix is defined 
as 


Rab{E) 


E 


(R) (R) 
ra<y Tab 

e]^'^ - E 


(23) 


{pi\HQQ\y) = H^, . (16) 

The sum extends over all open channels. The real shift 
function is defined by a principal-value integral. It is 
commonly assumed that the matrix elements Wafi change 
slowly with energy on a scale defined by the mean level 
spacing of the resonances. Then « 0. We use that 


This is in close analogy to Eq. (22), except that a fac¬ 


tor has been absorbed by each of the reduced width 
amplitudes. The resulting form of the scattering matrix 


IS 


Sab = 


+ 2*v^ ({1 - R{L - B)}-^R)^^ . 


(24) 
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Here Pa is the penetration factor in channel a, the ma¬ 
trices B and L are diagonal with elements Ba and La = 
Sa + iPa, and the real entities Sa — Ba play a role that 
is analogous to that of the shift function A in Eq. (151. 
The diagonal matrices L and P depend only on chan¬ 
nel radius and wave number ka [25]. The boundary 
condition parameter Ba is often taken as Ba = —la |26j . 
with la the orbital angular momentum of relative motion 
in channel a. In the i?-matrix approach the phases (pa 
are caused by elastic scattering on a hard sphere of ra¬ 
dius Ta while in the approach of Ref. m they are elastic 
potential scattering phase shifts. 


C. Implementation of stochasticity 


To fully define the scattering matrices in Sec. IIB 
we need to determine the resonance parameters. This 
is done by introducing statistical assumptions, using 
random-matrix theory |27j as a guiding principle. The ac¬ 
tual procedure is somewhat different for the three forms 


of the scattering matrix in Sec. IIB These are referred to 
as the S'-matrix approach, the it'-matrix approach, and 
the i?-matrix approach, respectively. 

The relevant random-matrix ensemble is the time- 
reversal-invariant Gaussian Orthogonal Ensemble 
(GOE). The elements Pl]SP^^ of the iV-dimensional 
GOE matrix are Gaussian-distributed real 

random variables with zero mean values and second 
moments given by 


\2 

= pSpaSup) ■ (25) 


(GOE) tr(GOE) 
tin 


H, 


Here and in what follows, the ensemble average is de¬ 
noted by an overbar. The parameter A is related to the 
average level spacing d at the center of the GOE spec¬ 
trum hy d = ttX/N. Universal properties are analytically 
derived m in the limit N —>■ c». These are: The eigen¬ 
values and the eigenvectors of statistically 

uncorrelated. The projections of the (real) eigenvectors 
on any fixed vector in Hilbert space have a Gaussian dis¬ 
tribution centered at zero. The eigenvalues obey Wigner- 
Dyson statistics. The degree to which these properties 
can be implemented depends on the approach used. 

In the ^'-matrix approach the Q-space Hamiltonian 


p[^i, of Eq. (|18|) is replaced by . That replacement 


provides the most direct implementation of random- 
matrix theory into scattering theory. The average cross 
section is worked out as an average over the GOE. For 
the specification of the parameters Wap one uses the in¬ 
variance of the GOE under orthogonal transformations 
in Hilbert space. That invariance implies that ensemble 
averages can depend on the lU’s only via the invariant 
forms J2p WapWpb- For the average S matrix to be diag¬ 
onal, the sums must be diagonal in the channel indices. 


and the only parameters left are the With 


(27) 


these determine the average S matrix elements Saa and 
the transmission coefficients Ta as 


Saa = 


1 - 

1 + ’ 


T = 

J- - 


4a:a 


(1 -f Xa)'^ 


(28) 


Equations (281 imply that the average strength Xa of the 


coupling of the GN resonance states to channel a is fixed 
by the average S matrix and, thus, determined by the 
shape-elastic input. In that sense, the GOE ensemble 
average of parameter-free. This is different from 

past calculations using a statistical R matrix or K ma¬ 
trix. 

In the itT-matrix and i?-matrix approaches, two as¬ 
sumptions are made: 

1. The partial width amplitudes yaa and the reduced 
width amplitudes yPa- both have a Gaussian dis¬ 
tribution with zero mean and a specified second 
moment ( 7 g). 

2. The eigenvalues Ea and Ea ' obey Wigner-Dyson 
statistics. 

If fully implemented, these assumptions correspond for 
iV —oo to the properties of the GOE listed below 


Eq. (25). 


In practical calculations, the implementation of these 
statistical assumptions causes difficulties. The S'-matrix 
approach lends itself to an analytical calculation of 
in the limit N ^ oo. The resulting expression is given in 


Eq. (35) below. However, the use of that expression was 


limited for a long time because of the difficulties in cal¬ 
culating reliably the ensuing threefold integral. A direct 
implementation would consist in drawing the elements 
from a Gaussian distribution, choosing a set of 


matrix elements Wap consistent with Eqs. (26) and (28), 


and inverting the resulting matrix Dp^. For N 3> 1 that 
is quite cumbersome and, to the best of our knowledge, 
has not been done before. We report on such a calcula¬ 
tion below. 

For the AT-matrix and i?-matrix approaches, it is 
straightforward to draw the partial width amplitudes or 
the reduced width amplitudes from a Gaussian distribu¬ 
tion. To meet postulate 2, the eigenvalues should be de¬ 
termined by diagonalization of the GOE matrix 77 (^ 10 ^) 
for N ^ 1. This is cumbersome, and a simplified version 
sometimes replaces postulate 2. The Wigner surmise for 
the distribution P(s) of spacings s of neighboring eigen¬ 
values reads 


Pw[s) = -sexp I- — 


(29) 


^ WapWpb = SabNvl , 


(26) 


with s the actual spacing in units of d. Spacings of neigh¬ 
boring eigenvalues are drawn at random from Pw{s) and 
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are used to construct the spectrum. Higher correlations 
between eigenvalue spacings are thereby neglected. In 
particular, the stiffness of the GOE spectrum (a central 
property) is not taken into account. 

With this input, the energy-averaged cross section 
(|(5a6 — Sab\‘^) can be calculated. It is often assumed that 
the energy average can be replaced by an ensemble av¬ 
erage over the joint distribution of level energies and de¬ 
cay amplitudes. The ensemble average can be readily 
obtained even in the limit of isolated resonances. The 
energy-averaged S matrix is more simply obtained using 
a Lorentzian average of width I and given by 


{S{E)) = S{E + il) , (30) 


and S{E + il) is obtained by replacing K{E) in Eq. (22) 
or R{E) in Eq. (23) by K{E + il) or R{E + il), re¬ 
spectively. That yields the transmission coefficients in 


Eq. ([5|. 

Monte Carlo calculations based on this approach ^5] 
have been used to define heuristic parametrizations of 
the width fluctuation correction factor Wat- The re¬ 
sults of Hofmann, Richert, Tepel, and Weidenmiiller [8l 
[U (HRTW) are based on the K matrix, those of 
Moldauer m on the R matrix. The resulting fit for¬ 
mulas for the WEC factor are collected in the Appendix. 



Incident Neutron Energy [MeV] 



FIG. 1. (Color Online) An example of an elastic scattering 
cross section generated with the statistical R-matrix method 
(bottom panel), compared with the real cross section in 
ENDF/B-VII.l (top panel), for a neutron-induced reaction 
on ®®Fe in the 1-2 MeV energy range. 


III. MONTE-CARLO SIMULATIONS 
A. R matrix and S matrix 

In the 1960s and 70s, the statistical i?-matrix approach 
used by Moldauer H [S] offered the only possibility to 
use random-matrix ideas in CN scattering. As an ex¬ 
ample for that method we show in Fig. the result of 
a new Monte-Carlo simulation of the elastic cross sec¬ 
tion for neutron scattering on ®®Fe (bottom panel). This 
is compared with the real cross section (upper panel) 
given in ENDF/B-VH.l [21]. In the simulation we put 
(7c) = 10 keV for the s wave, and 200 eV for the higher 
partial waves (p, d, and / waves). The radiative capture 
channel was ignored. Although the statistical i?-matrix 
calculation cannot reproduce the detailed structure, the 
comparison provides information on average properties 
and, thus, a useful link between the fluctuating cross sec¬ 
tion and the optical model calculations |5D] . 

However, the approach involves a number of parame¬ 
ters, such as the partial widths, the level density EH ED, 
the energy range of interest, and so forth. Strict imple¬ 
mentation of the more abstract GOE approach actually 
removes the need to define these parameters. This is 
most easily demonstrated for the case of the K matrix. 
Upon scaling the energies by the mean level spacing d so 
that E/d e, E^/d Ca, the quantities e and are 
dimensionless, and the spacing distribution of the Co- is 
given in terms of the universal dimensionless correlation 
functions of the GOE m- The expression for Pw{s) in 


Eq. (29) is an example. Applying the analogous scaling to 
the partial width amplitudes, 'jaa/d^^'^ —)■ jaa generates 
dimensionless uncorrelated Gaussian-distributed random 
variables The second moments of these quantities 
are determined by the average iS-matrix elements of the 
optical model. In other words, the scaling E/d —>■ e, 
E^/d ea, laa/d^^"^ —>■ 7aCT maps the GN scattering 
problem onto a GOE scattering problem where the only 
input parameters are the number of open channels and 
the elements (Saa) of the average scattering matrix in 
each channel. That scattering problem describes univer¬ 
sal chaotic scattering. 

For the R-matrix simulation, Eqs. 0 and ( [I^ with 
replaced by as 

W,b , (31) 

flU 

Df,, = E5^, - , (32) 


serve as starting point. The matrix elements obey 
Eq. (26). The simulation generates an ensemble of S ma¬ 


trices by generating a single set of matrix elements {W^a} 
combined with a number of realizations of By 

construction, all S matrices in the ensemble have for 
iV —>■ oo the same mean values. The matrix elements 
Wfia are determined as follows. Given the elements G^a 
of a coupling strength matrix G of dimension N x A, 
where A is the number of channels, diagonalization of 
the real symmetric matrix G^G in channel space with 
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an orthogonal matrix O yields 

O-^G^GO = di&g{vl) 
W = GO . 


(33) 

(34) 


This procedure guarantees that Eqs. (26) are satisfied. 
The eigenvalues and d = nX/N define Xa = n'^vlfd 
and these, in turn, the transmission coefficients Ta via 
Eq. (28). Figure shows examples of calculated elastic 
scattering cross sections for A = 2, TV = 20,100, and 
for three different transmission coefficients Ta = 0.1, 0.5, 
and 0.99. To show how the cross section evolves as the 
transmission coefficient increases, we fixed the random 


number sequence so that the eigenvalues of are 

the same for the three r„ cases. 


B. Ensemble average 


The ensemble average of Eqs. © and ( [3^ can be eval¬ 
uated numerically either by employing the MC technique 
where the elements from a Gaussian 

distribution, or by calculating the three-fold integral of 
Verbaarschot, Weidenmiiller, and Zirnbauer m 


_ -| poo poo pi 

= -1 dX,J^ MA,Ai,A2) 


X e 


0 ^0 

—ir(Ai +A 2 + 2 A) 


i-t;a 


V(1 + T,Ai)(1 + T,A2) 


J(A, Ai, A 2 ) , 


A(1 — A)|Ai — A 2 I 


where 

/r(A, Ai, A 2 ) = 

VAiti -h AijV'^2fi + 

J(A, A„ A,) = + j^) 

A2(1 -I- A 2 ) 


vA7(T^i-Aiy-\/A2Xi~-fA^(A -I- Ai)2(A + A2)^ 

-+ 

l+Ta) 

Ai(l -I- Ai) 


{dacdbd T daddbc)TaTb 


r = -{E2-Ei) . 


2X{1 + X) 


(l+TaXlKl+TbXl) il+TaX2)il+TbX2) {l-TaX){l-TbX) 


(35) 

(36) 

,(37) 

(38) 


The triple-integral in Eq. (35) can be evaluated numer¬ 
ically by introducing new integration variables |33j that 
avoid singularities in the integrand, and by the Gauss- 
Legendre quadrature with the order high enough to ob¬ 
tain convergence m- I n pr actical applications we need 
SabS*f^ only, so that Eq. (351 can be reduced to a slightly 
simpler form [2T] as SaaSaa = 1 — Ta- This is not the 
case if we have off-diagonal elements in (S'), or different 
energy arguments so that r ^ 0. 

A benefit of MC is that we are able to explore a 
larger parameter space, while Eq. (35) holds in the limit 
N ^ 00 . In Eq. (35) we have replaced the energy average 


(|(5a6-S^°°^^P) by the ensemble average \6ab - 
The difference between the two averages is discussed 
later. Hereafter we always calculate the ensemble average 
unless stated explicitly otherwise. The average is evalu¬ 
ated at the center of the GOE eigenvalue distribution, 
E — 0. As shown in Fig. the calculated cross section 
near E = 0 for a single realization of jispiayg 

chaotic fluctuations. The number of MC realizations 
needed to obtain a meaningful average varies from 10,000 
to a million, depending on convergence. The criterion 
used was that the deviation of the average S'-matrix from 


its input value was sufficiently small, lA^abl < 10“®. 

Figurej^shows the probability distribution of the elas¬ 
tic scattering cross section at if = 0 for A = 100, A = 2, 
and three different Ta = Tb values of 0.1, 0.5, and 0.99. 
The MC ensemble average values are indicated by the 
location of the arrows, e.g., in the case of Ta = 0.99, the 
average is 1.47. To compare these averages with predic¬ 
tions of the statistical model, we have to subtract the 
direct part (1 — JftS'aa)^ from the elastic channel. That 
gives the average fluctuating part of 0.660. The 

Hauser-Feshbach cross section is 

= 0.495 , (39) 

giving in that case a 25% smaller elastic cross section. 
When the GOE triple-integral of Eq. ( [3^ is performed for 
the given transmission coefficients, the simulated cross 
sections are recovered. In Table [l] we compare the MC 
results with other statistical models — KKM [5] , HRTW 
mm, Moldauer [T3], Ernebjerg and Herman [52], and 
Kawano and Talou [23]. In general, all the statistical 
models predict the average reasonably well when {T)/d 
is large. More comparisons of the MC generated cross 
sections with these statistical models can be found in 
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FIG. 2. (Color Online) Generated elastic |1 — S'aap and inelastic scattering cross sections with the GOE 5'-matrix for 

three different transmission coefficients Ta of 0.1, 0.5, and 0.99. The solid curves are the elastic, and the dotted curves are the 
inelastic cross sections. The left column is for A'’ = 20 and the right is for the N = 100 case. 


Ref. p5] . 

One may argue that the agreement between the GOE 
triple-integral and the MC simulation is obvious because 
the triple-integral is an analytical form of the ensemble 
average for Eqs. (31) and (32) in the limit of iV —?> oo. 
We have, therefore, studied the fV-dependence of the cal¬ 
culated averages. Starting with N = 100, we reduce the 
number of resonances and compare the ensemble aver¬ 
age of the elastic cross section with the triple-integral 
results. Surprisingly the triple-integral still gives very 
accurate average values even if TV = 3. Averaging over 
a few resonances is certainly an extreme case, and is not 
realistic. 


IV. VALIDATION OF STATISTICAL MODELS 
A. Energy average versus ensemble average 


ble average can be performed numerically using the MC 
simulations for finite N^ and (c) the average is taken over 
energy and calculated for a single realization of the en¬ 
semble. Method (c) is the only way to perform averages 
over actual data. Such averages define the optical model. 
Obviously it is highly important to know whether (and 
if so, when) these averages agree. 

Let w{Eo, E, I) be the weight function centered at en¬ 
ergy Eq with width I used to define the average over 
energy E. In what follows w{Eq, E, I) is taken to be a 
Lorentzian. Our aim is to know under which circum¬ 
stances the equality 




w{Eo, E, I)S{E)dE = S{Eo) 


(40) 


holds. Since there is no analytical way to investigate that 
relation, we ask when the weaker condition 


There are three ways to calculate averages: (a) the 
ensemble average can be performed analytically in the 
limit N —> oo, which is given in Eq. (35), (b) the ensem- 


|(5)-S'| =0 (41) 

is fulfilled [Tl]. It is straightforward to show that Eq. (41) 


















































































































TABLE I. Comparison of numerical average IS'aaP for some cases of Ta = Tb = 0.1, 0.5, and 0.99, with the statistical models 
— Hauser-Feshbach [^, KKM [^, HRTW [11], Moldauer [T3|, GOE |15| . Ernebjerg-Herman [22], and Kawano-Talou |23| . 


Ta 

0.1 

0.5 

0.99 

Elastic Inelastic Elastic Inelastic Elastic Inelastic 

MC simulation 

0.0733 

0.0261 

0.351 

0.149 

0.660 

0.330 

Hauser-Feshbach 

0.0500 

0.0500 

0.250 

0.250 

0.495 

0.495 

KKM 

0.0662 

0.0332 

0.333 

0.167 

0.660 

0.330 

HRTW 

0.0737 

0.0257 

0.352 

0.147 

0.661 

0.330 

Moldauer 

0.0734 

0.0260 

0.349 

0.150 

0.665 

0.325 

GOE 

0.0734 

0.0260 

0.351 

0.148 

0.661 

0.330 

Ernebjerg-Herman 

0.0742 

0.0252 

0.366 

0.134 

0.681 

0.310 

Kawano-Talou 

0.0735 

0.0259 

0.351 

0.148 

0.661 

0.330 



We assume that the level spacing d = nX/N is indepen¬ 
dent of energy. We limit ourselves to the case where all 
transmission coefficients are equal and given by T^. We 
perform the energy averages using Lorentzians centered 
at zero, 

= (44) 

with width I specified in units of d/ir. We define 

/ OO pOO 

/ w{Ei,0,I)w{E2,0,I) 

-OO j —OO 

X R{E2 - Ei;Ta, A) dE^dE^ , (45) 

where 


RiE2-Ei-TaA) = — - - -(46) 

V 2 1, a, ; ^ ^ 


FIG. 3. (Color Online) Distribution of elastic scattering cross 
section at iJ = 0 for many GOE 5'-matrix realizations. Three 
cases, Ta = Tb = 0.1, 0.5, and 0.99 are shown. The arrows 
show the actual average values for each distribution. 


is equivalent to 


r+oo 




dEiw{Eo,Ei,I) 


dE2w{Eo,E2,I) 


S^{Ei)S^*{E2)=0 , 


(42) 


where the two-point function S^{Ei)S^*{E 2 ) is given by 
Eq. (|3^. 


The average two-point function S^{Ei)S^*{E 2 ) in¬ 
volves two S matrices at energies Ei and E 2 . Because of 
the very weak energy dependence of the average S ma¬ 
trix we approximate S{Ei) ~ S{E 2 ) and evaluate both 
at E = 0. Then the energies Ei and E 2 in the two-point 
function appear only in the oscillating term 


exp{-fr(Ai-I-A 2 + 2A)} , 


^{E2-E,). (43) 


We use the real part only because integration over the 
imaginary part in Eq. ( |4^ yields zero. 

Our results for the elastic channel are displayed in 
Figs. |^to[^ Figure 1^ shows the function L of Eq. (45) 
versus I for A = 10 and for values of Ta ranging from 
0.1 to 0.9. As expected, L decreases as / increases so 
that ensemble average and energy average agree when 
the Lorentzian width I is sufficiently large. To get the 
same accuracy larger values of Ta require larger widths 
I. The dependence of L on channel number A is shown 
versus I in Fig. for Ta = 0.5 (logarithmic scale) and 
in Fig. I^for Ta = 0.99 (linear scale). Larger values of A 
require larger values of /, the slowest decrease occurring 
for the strong-absorption case where Ta is close to unity. 
For the strong-absorption case Ta — 0.99, Fig. shows 
the values of I versus channel number for which L = 0.1. 
The result is a clear linear dependence 


/(L = 0.1) ~ 2.2A-b 1.9 . (47) 

In the Ericson regime J2a 1 or, for equal trans¬ 

mission coefficients in all channels, T^A ^ 1, the auto¬ 
correlation function is known analytically. The real part 


r = 
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Lorentzian Width, I 


FIG. 4. (Color Online) L{Ta, A, I) as defined in Eq. (451 
versus the Lorentzian width I for A = 10 channels and for 
different values of the transmission coefficient Ta- From the 
lowest to the highest cnrve Ta changes from 0.1 to 0.9 in steps 
of width 0.1. 


FIG. 5. (Color Online) L(0.5, A, 7) as dehned in Eq. ( |45[ ) ver¬ 
sus the the Lorentzian width I for Ta = 0.5 and for different 
channel nnmbers A. From the lowest to the highest curve the 
values of A are 2, 5, 10, 20, 30, 40, and 50. 


is a Lorentzian with denominator where the to¬ 
tal width is given by F = (d/27r) large I the 

function L falls off with (2/)“^. We have L = 0.1 for 
I « 5F. 

The rate of decrease of L versus I depends on Ta and 
A. Using our results we can nevertheless draw some gen¬ 
eral conclusions concerning neutron-induced reactions at 
low energy. In the domain of isolated resonances the 
number of channels is effectively small (7 channels are 
numerous but extremely weak individually). Flere the L- 
function becomes ~ 0.1 or less when I is larger than 10 
or so. A value of 7 = 100 corresponds to lOOd/Tr ~ 30d. 
Hence the L-function will be sufficiently small when the 
energy-averaging interval is one or two orders of mag¬ 
nitude larger than the average resonance spacing d. In 
the Ericson regime that same statement applies with d 
replaced by F, the average total resonance width. 


B. Asymptotic value at strong-absorption limit 

1. Elastic enhancement factor in Ericson limit 



0 5 10 15 20 25 30 35 40 

Lorentzian Width, I 


FIG. 6. (Color Online) L(0.99,A, 7) as defined in Eq. (451 
versus the Lorentzian width 7 for Ta = 0.99 and for different 
channel numbers A. From the lowest to the highest curve the 
values are A = 2, 3, 4, 5, 6, 8, 10, 12, and 15. The horizontal 
line shows L = 0.1. 


In the strong-absorption or Ericson limit 
Eq. (35) yields Wa = 2 for the elastic enhancement factor 
or, equivalently, Va = 2 for the channel degree-of-freedom 
|33j . Explicitly we have 


(CTah) = 


(1 + 5ab)Tan 

EcTc 


(48) 


The dots indicate terms of order (X^c^c) ^ or higher. 
The term of leading order is the Hauser-Feshbach result 


with an elastic enhancement factor of two. Most sta¬ 
tistical models agree with that result. An exception is 
the model by Moldauer, which has an asymptotic value 
of izq = 1.78. Although Moldauer’s heuristic method 
to obtain Eq (A.7) in the Appendix is somewhat simi¬ 
lar to the MC technique we adopt here, there is a no¬ 
table difference between the two approaches. In the MC 
approach we perform the ensemble average over the ele¬ 
ments of the Hamiltonian 77^^^^^ Moldauer’s statistical 
77-matrix model has two independent inputs: the decay 














10 



FIG. 7. (Color Online) The symbols show for which values of 
the Lorentzian width I and channel number A the function 
L{0.99, A, 1) attains the value 0.1. The line is a least-squares 
fit to the symbols; I{L — 0.1) ~ 2.2A + 1.9. 


widths drawn from the Porter-Thomas distribution, and 
the level spacing sampled from the Wigner distribution 
in Eq. (29). 



FIG. 8. (Color Online) Distribution of decay amplitudes 'jaa, 
when the GOE S matrix is written in the AT-matrix form. The 
histograms correspond to Ta = 0.1, 0.5, and 0.99, respectively. 
For Ta = 0.99, we compare the Gaussian distribution with the 
one obtained from the standard deviation in Eq. (59). 


2. Decay amplitude distribution 


Our aim is to reproduce Moldauer’s lower asymptotic 
value by modifying the MC sampling method. Before 
doing that, we show the distribution of the width ampli- 
tudes v^Tao- when we rewrite our stochastic S'-matrix of 
Eq. (31) in an equivalent form [niiiH] 


WaaWab 


(7 

V 


(49) 

(50) 

(51) 


where is the eigenvalue of form 

the width amplitudes Waa = are uncorrelated 

Gaussian-distributed random variables with zero mean 
values and the standard deviation. We produced the dis¬ 
tributions of W for the case N = 100, A = 2, and three 
values of Ta = 0.1, 0.5, and 0.99. 

The width distributions are shown in Eig.j^for the elas¬ 
tic channel. Because we used the same transmission for 
both channels, the distribution for the elastic and inelas¬ 
tic channels are identical. Figure shows the standard 
deviation Ua for each Gaussian for various Ta- 

The second moment of Gaussian distribution is given 
by m 


which is shown by the two dashed curves in Fig. The 
sign ambiguity in Eq. ( [^ is caused by the fact that there 
are two values of Saa with opposite signs that yield the 
same value of Ta = 1 — |S'aaP- 


3. Emulating Moldauer’s calculation 


Moldauer’s K matrix (see Section IIB) can be written 


as 


K^,{E) = 5ab^Kl + Y. 

<7 


E-E, ■ 


(53) 


The elements of the elastic background matrix and 
the variances of the amplitudes w are determined by the 
energy-averaged S matrix. For TT® we have 


1 + 5(goE)(^ + j7) 


(54) 


showing that is determined by the transmission co¬ 
efficient Ta. Since QS{E + il) ^ 0, we may omit the 
background term When we view the TT-matrix as 

an i?-matrix, QK(E -|- il) is the pole strength 2'K^'^/d, 
therefore the second moment for the distribution of the 
widths Waa reads 


a 


2 

a 


d 


Ta 


(52) 


o-a = 27r72 = - \QKl 


7T 2 -Ta± v/r^ 


TT 


(55) 
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Tranmission Coefficient 


FIG. 9. (Color Online) Standard deviation of GOE decay am¬ 
plitude distribution as a function of transmission coefficient. 
The dashed curves are Eq. (521 and the solid curve is Eq. (591. 


The elastic enhancement factor Wa can be defined only 
when all channels are identical, Ta = Th = ... = T\. 
That is the case we address. 

The calculation of the ensemble average of Eq. (53) 


proceeds as follows. First we generate S of Eq. (31), and 


convert it into K via Eq. (21). As Moldauer performed in 
Ref. [13], we use K{E + ir)&nd Eq. (55) to determine the 
average widths of the decay amplitudes. The latter are 
then sampled from Gaussians with widths Ua-, indepen¬ 
dently of the GOE eigenvalues. The Lorentzian average 
width I is taken to be 0.2 A. We extract the elastic en¬ 
hancement factors and compare with the standard GOE 
simulation that is described in Sec. ttUBi 

The elastic enhancement factor Wa is calculated as 


Waa = 


rHF 


rp2 

HF _ 
aa 


a 


Wa = 


(A - l)Wa 
' A - Waa 


Wa-1 


(56) 

(57) 


We calculate K{E + il) for each realization of the GOE 
^-matrix. Therefore, the ensemble average of Eq. (53) 


converges slowly. In addition, simulations for very large 
values of or A are not feasible in general. We chose 
N = 200, A = 5, 10, 20, and 30. The transmis¬ 
sion coefficients are 0.25 and 0.75. These combinations 
roughly cover Moldauer’s numerical study of the strong- 
absorption cases. 

The values of Va versus obtained in that way 

are compared with the GOE result in Fig. [T^ The sym¬ 
bols in the upper panel show the results of the standard 
GOE simulation, those in the lower panel the results of 


the simulation described in the previous paragraph. The 
curves in the upper panel represent Eqs. ( |A.19 ), those in 
the lower panel represent Eq. (A.7), both for the cases 
Ta = 0.25 and 0.75. These equations are meant to ap¬ 
proximate Va for given values of the transmission coeffi¬ 
cients. The results of the GOE simulation are well rep¬ 
resented by Eq. (A. 19) which has the asymptotic value 
of 2 in the strong-absorption limit. The MG simulation 
that uses Eq. (55) tends to give lower Va values, similar 
to Moldauer’s findings. 

A plausible explanation of this discrepancy relates to 
the determination of the decay amplitude via Eq. (55). 
Since 


= Y^a ’ 

the widths in Moldauer’s approach have a second moment 
given by 


d Ta 


TT 2 - T„ 


(59) 


which is shown in Fig. by the solid curve. Comparison 
with Eq. (52) shows that this is correct only for small 
values of T^. Discrepancies arise for Tq « 1. In Fig 
we compare for Ta = 0.99 the distribution of widths us¬ 
ing for the second moment the correct expression (|52[ ) 
with the one obtained from Moldauer’s equation (59). 


(We do not show the Ta = O.I and 0.5 cases because 
they perfectly overlap with the exact values). We note 
that Moldauer’s approach gives a slightly narrower dis¬ 
tribution. We suspect that this is the root of Moldauer’s 
incorrect asymptotic value for Va = 1.78. 


4- Asymptotic expansion 


The next-to-leading-order term of Eq. (48) is given by 


an asymptotic expansion of Eq. (35) in inverse powers of 
EcTc 


], which is also given in Appendix. This is 
shown by the dashed curves in Fig. 10 (a). The asymp¬ 


totic expansion approximates the GOE triple-integral 
very well, when Yc ^ This might be practically 
useful in the strong-absorption limit, in particular when 
the number of open channels is so large that calculation 
of the GOE triple-integral becomes extremely difficult. 


C. Very weak entrance channel 

An extreme case where all the statistical models fail is 
reported in Ref. |23|. When there are few open channels 
with either very small or very large transmission coeffi¬ 
cients, none of the width fluctuation models reproduces 
the GOE results. That was also discussed by Moldauer 
[35] as the total width fluctuation, and his numerical 
study shows a strong enhancement in the elastic chan¬ 
nel. We performed the GOE simulation for the case of 
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FIG. 11. (Color Online) Ratio of the elastic cross section 
to the Hauser-Feshbach prediction as a function of the ra¬ 
tio Ta/Th of transmission coefficients. The symbols are the 
MC simulation results. The curves are predictions by various 
statistical models. 


FIG. 10. (Color Online) Channel degree-of-freedom values Va 
as functions of Ta- The symbols are the Monte-Carlo sim¬ 
ulation results, see text. The solid curves in the top panel are 
from Eq. ( A.17[ ), and the dotted curves are from Eq. (A. 101 
for Ta = 0.25 and 0.75. The curves in the bottom panel show 
Moldauer’s systematics given by Eq. (A.71 for the same set of 
Ta. 


= 100, A = 2 and Ta/Tb ^ 1. The calculated width 
fluctuation correction factor Waa, which is the ratio of 
the elastic channel cross section to the Hauser-Feshbach 
cross section, is shown in Fig. El Since the GOE triple¬ 
integral is correct for all values of A and Ta, the MC 
simulation perfectly agrees with GOE, except some devi¬ 
ation seen at very small Ta/Tf, values, due to numerical 
instability. 

Few-channel cases with very different values of the 
transmission coefficients are very special and hard to re¬ 
alize in practice. A photo-induced reaction that creates 
a compound nucleus just above neutron threshold could 
be a case in point. However, since almost all incoming 
flux goes to the neutron channel and to the other gamma 
channels, the photon compound elastic cross section is 
tiny even if it is enhanced by a factor of 50. That is why 
it might be difficult to confirm the strong enhancement 
in the elastic channel experimentally. 


V. DIRECT REACTIONS 
A. Engelbrecht-Weidenmuller transformation 

So far it was assumed that the average S matrix is di¬ 
agonal. That assumption fails when some channels are 
strongly coupled. In practice that happens, for instance, 
when collective states in the target nucleus are excited 
by an incident nucleon (a direct reaction). In such cases, 
the average S matrix is not diagonal. The unitarity of 
the scattering matrix imposes strong constraints on the 
scattering amplitudes. As a consequence, directly cou¬ 
pled channels cause correlations between the resonance 
amplitudes in those channels. That is why the calcu¬ 
lation of the average compound-nucleus cross section in 
the presence of direct reactions has been a long-standing 
problem. 

When (S) is not diagonal, the definition of the trans¬ 
mission coefficients T must be generalized. That is done 
using Satchler’s transmission matrix 

Pab = Sab-Y.(^--)(^bc) ■ ( 60 ) 

c 

In the strong-absorption limit, Kawai, Kerman and 
McVoy (KKM) [B] expressed the compound-nucleus cross 
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section in terms of the matrix P (see Eqs. (A.8) and 
( |A.9 1). Actual calculations using KKM including the 
direct channels are, unfortunately, very limited, e.g. 
Refs. [57] and |55] . 

In practical calculations, an often-used approximate 
way to include the direct reaction in the statistical model 
consists in redefining the transmission coefficients so as 
to take account of some direct reaction contribution. 




(61) 


A closed form of the average cross section based on the 
GOE triple-integral formula that takes the EW trans¬ 
formation into account, was derived by Nishioka, Wei- 
denmiiller, and Yoshida [43] . However, the computation 
might be impractical. We follow the EW transformation 
step-by-step from Eq. (601 to Eq (65). The result allows 


us to estimate uncertainties due to the approximation 


Eq. (61). 


B. Ensemble average using EW transformation 


The sum of the modified transmission coefficients T' 
equals Tr(P). Therefore, it is reasonable to expect that 
GOE cross-section calculations using the modified trans¬ 
mission coefficients as input parameters as done in 
Ref. [39] may not be far off the mark. In comparison with 
the exact approach introduced below, the method greatly 
simplifies the calculations. However, a quantitative vali¬ 
dation of the simplification (611 and an understanding of 


its limitations are still needed. 

The following rigorous treatment of the direct reaction 
was proposed by Engelbrecht and Weidenmiiller (EW) 
goj. Since P is hermitian, P can be diagonalized by a 
unitary matrix 


{UPU^)ab = SabPa , 0 < Pa < 


< 1 


The transformation U also diagonalizes the average scat¬ 
tering matrix, 

{S) = U{S)U^ with {S)ab = Sab{S)aa ■ (63) 


To implement direct reactions, one may use, for in¬ 
stance, the pole expansion of the S matrix. We find it 
simpler to employ the AT-matrix as in Eq. (22). We allow 
for a direct background by writing 


K,b{E) = + 


WaaWab 

^ E-E„ 


( 66 ) 


where the elements of the background matrix serve 
as parameters. When K is real and symmetric, S is au¬ 
tomatically unitary. 

We consider a case with direct coupling between two 
channels only. The background matrix is 


^( 0 ) 


kaa kab 0 . . . ^ 

^ab kbb 0 ••• 

0 0 0 ••• 


V ; : ; y 


(67) 


In the diagonal basis of P, the transmission coefficients 
are given by 


Pa = 1 - \{Saa)\^ ■ (64) 

In that basis, the decay amplitudes in different chan¬ 
nels are statistically uncorrelated, and the calculation of 
SpqS*g proceeds as described above for the case without 
direct reactions, with pa as input parameters. The result 
must be transformed back to the physical channels. That 
gives |S] 


= E u;,u;,UraUsbSpqS*s • (65) 

pqrs 

Moldauer demonstrated the impact of the EW trans¬ 
formation numerically |5]. He argued that the flux 
into the strongly coupled inelastic channels is enhanced. 
Gapote et al. [IT] demonstrated that enhancement by 
applying the coupled-channels code EGIS [12 to neutron 
scattering off Although EGIS is capable of perform¬ 
ing the EW transformation, it has some approximations 
and limited functionality, particularly for calculating the 
neutron radiative capture and fission channels. The EW 
approach uses only the average S matrix as input and 
facilitates showing how direct reactions impact on the 
compound nucleus. 


Eor the sake of simplicity, we take kaa = kab 
where fcp is real. The average S matrix is 


- _ 1 - + TTjWaWb) 

1 -h iATW - TT{WaWb) ' 


kbh — ko , 


( 68 ) 


The amplitudes Waa are zero-centered Gaussian- 
distributed random variables, uncorrelated for a ^ b. 
The parameters then are N, A,Ta, k^. For simplicity we 
use the same for all channels. 

The cross sections are calculated in the following three 
ways. 

• For each value of ko, the MG method is used to 
generate 100,000 realizations of S. The average 
cross section is obtained directly as the average of 
ISab — Sabl"^ over that ensemble. Figure [T^ shows the 
cross sections for N = 100, A = 2, Ta = 0.8 = Tt, 
and for fcg varying from 0 to 2 obtained in that way. 
The top panel shows the elastic scattering cross sec¬ 
tion |1 —S'aaP, the bottom panel shows the inelastic 
scattering cross section l^afep- 


• The average of S over the ensemble of 100,000 re¬ 
alizations is used to calculate the modified trans¬ 


mission coefficients of Eq. (611. These are used 
in the GOE triple-integral to calculate the width 
fluctuation correction. 
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• S as obtained in the previous step is diagonalized 
using the EW transformation. The eigenvalues Pa 
are used in the GOE triple-integral. The result 

SpqS*g is back-transformed to |S'a6p. 

We analyze our results in terms of the usual “optical 
model” cross sections 


CTT = 2(1 - , 

CrSE = |1 — 5'aaP , 
O’DI = . 


(69) 

(70) 

(71) 


Here ctt, ctse and (Tdi stand for the total, the shape 
elastic, and the direct inelastic cross section, respec¬ 
tively. The reaction cross section and the compound 
formation cross section are defined as ctr = ctx ~ o'se 
and CTcN = ctr, — (Tdi, respectively. All these cross sec¬ 
tions are given by the coupled-channels optical model, 
while the compound elastic (uce) and compound inelas¬ 
tic (cTci) cross sections require statistical-model calcula¬ 
tions. We do not use a coupled-channels optical model 
in the present context but are able to calculate all these 
cross sections directly from the MC simulation. The pa¬ 
rameter fco controls the strength of aui up to a limit 
dehned by unitarity — since (Tt and (tse are connected 
by Saa, o’R is constrained even if fco is very large. 

Figure [T3| shows how the compound-inelastic scattering 
cross section changes with the strength of the direct reac¬ 
tion. We plot the ratio of aci to the reaction cross section 
ctr as a function of the ratio ctdi/o’r. The upper panel 
is for Ta — 0.5 and the lower panel is for = 0.9. In 
each panel we show two cases, A = 2 and 10. The results 
from the EW transformation agree perfectly with the MC 
simulations, confirming that the EW transformation with 
the GOE triple-integral yields the correct average cross 
section when there are strongly coupled channels. 

When the background AT-matrix is parametrized as in 
Eq. (67), CTci approaches the unitarity limit for very large 


fco. At this limit, we have (tce — o’ci. The elastic en¬ 
hancement disappears when a direct channel becomes 
very strong. Since we employed the same transmission 
coefficients for all channels, the compound elastic and in¬ 
elastic scattering cross sections are equal in that limit and 
given by ctcn/A. Use of the modified transmission coeffi¬ 
cients overestimates ctce and underestimates ctci. The 
discrepancy increases with increasing (Tdi. 

The EW transformation is definitely required to cal¬ 
culate the correct compound cross sections when A is 
small and (Tdi/(Tr is larger than about 5%. A case in 
point might be a reaction induced by neutrons of sev¬ 
eral 100 keV impinging on an actinide. Several levels 
of the ground-state rotational band will be excited by 
the direct inelastic scattering process. A simple coupled- 
channels calculation for the 300-keV neutron-induced re¬ 
action on 238 u gives ctdi/(Tr of about 0.1. Therefore the 
approximate method that uses the modified transmission 
coefficients is expected to result in an underestimate 
of (TCI. 



Energy [in unit of lambda] 



Energy [in unit of lambda] 


FIG. 12. (Color Online) Simulated elastic (top panel) and 
inelastic (bottom panel) scattering cross sections as functions 
of the backgronnd parameter fco. 


VI. CONCLUSION 


We have investigated the statistical properties of the 
scattering matrix containing a GOE Hamiltonian in the 
propagator. That S matrix describes general chaotic 
scattering and applies to compound-nuclear reactions at 
low incident energies (below the precompound regime). 
We have compared results for average cross sections ob¬ 
tained from Monte-Garlo (MC) simulations with those 
from the GOE triple integral and from statistical mod¬ 
els. The latter give heuristic accounts of the width fluc¬ 
tuation correction. In the GOE approach, the results 
depend on few parameters: the number N of resonances, 
the number A of open channels, and the average S ma¬ 
trix elements. Without direct reactions, the average S 
matrix is diagonal, and the relevant parameters are the 
transmission coefficients in the channels. When the chan¬ 
nels are strongly coupled and the average S matrix is not 
diagonal, the number of parameters is correspondingly 
increased. Our simulations indicate the range of valid¬ 
ity of the heuristic models and have led to the following 
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Direct Inelastic C.S. / Reaction C.S. 


FIG. 13. (Color Online) Ratio ctci/ctr of the compound in¬ 
elastic scattering cross section to the reaction cross section as 
a function of the ratio ctdi/ur of direct reaction cross section 
to the reaction cross section. The symbols show the ensem¬ 
ble average of the MC simulation, the dotted lines are cross 
sections calculated with the modified transmission coefficients 
Ta, and the solid lines are the result of the EW transforma¬ 
tion. The top panel is for Ta = 0.5, and the bottom panel is 
for Ta = 0.9. 


erage of MC-generated cross sections coincides with 
the result of the GOE triple-integral formula (351. 
Although that formula is derived in the limit of a 
large number of resonances, it gives the correct av¬ 
erage even if the number of resonances is small. 


• Energy average and ensemble average agree reason¬ 
ably well (i) for isolated resonances when the width 
of the Lorentzian averaging function is one or two 
orders of magnitude larger than the average reso¬ 
nance spacing and (ii) in the Ericson regime when 
the width of the Lorentzian averaging function is 
one or two orders of magnitude larger than the av¬ 
erage total width of the resonances. 

• In the strong-absorption limit (Ericson regime) 

where 1, the channel degree-of-freedom 

Va is 2, different from Moldauer’s asymptotic value 
of 1.78. 


• In extreme cases where a few open channels (includ¬ 
ing the incident channel) have very small transmis¬ 
sion coefficients and a few others have transmission 
coefficients close to unity, the elastic channel is sig¬ 
nificantly enhanced. Most of the standard statisti¬ 
cal models cannot predict that enhancement. The 
GOE triple integral is the only way to produce the 
correct average cross section. 

• Direct reactions (for instance, the excitation of 
states of a rotational band due to inelastic scat¬ 
tering) cause the average S matrix S to acquire 
large off-diagonal elements. Using the Engelbrecht- 
Weidenmiiller (EW) transformation we have diago¬ 
nalized S and evaluated the GOE triple integral in 
the diagonal channel basis. The results agree with 
the MG simulations. We find that the direct reac¬ 
tion increases the inelastic cross sections while the 
elastic cross section is reduced. 
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Appendix: Statistical models 
1. HRTW 


In the HRTW approach [SlE], an elastic enhancement 
factor Wa is expressed by the channel transmission coeffi¬ 
cient Ta, and all the channel cross sections are calculated 
from an effective transmission coefficient Va 

{<Jab) = {1 + - 1)} , (A.l) 

Z^C 


where VcS are determined from the unitarity of S'-matrix, 
in another word, the flux conservation. The values of 
Wa were derived from the statistical iC-matrix analysis. 
There are two sets of Wa parameterization, namely in the 
original paper of Ref. [8] , and the updated parameters in 
Ref. m- We refer to the updated parameters as HRTW, 
which reads 


Wa = l + 


F = 4 


T 

f 



(A.2) 

(A.3) 


where T is the average value of Ta, and T is the sum of 
Ta for the all open channels T = J2c 


2. Moldauer 


The Gaussian distribution of yields the Porter- 
Thomas distribution of j'^a when there is only one chan¬ 
nel. More generally, the distribution of will be the 
distribution with the channel degree-of-freedom i/q. 
In these circumstances, the width fluctuation correction 
factor can be evaluated numerically as IIISIES] 


Wab = (1 + 


2Sab 

I'a 


io Fa{t)Fb{t)WFk{tY^/^'^ ■ ^ 
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F,{t) = l + ~t. (A.5) 

The integration can be performed easily by changing the 
variable t into z as 


z dt 1 

1-z’ di ^ (l-z)2 ’ 


(A.6) 


where z —>■ 0 for t = 0, and z —>■ 1 for t —>■ oo. 

In contrast to HRTW, Moldauer’s prescription gives 
the width fluctuation correction factor that ensures the 
unitarity for all the channels when the channel degree- 
of-freedom Va is provided. Moldauer obtained Va as a 
function of each channel transmission coefficient and 
the sum of them T = with the MC simulation, 

which reads [T3] 


Va = 1.78 + - 0.78) exp(-0.228T) . (A.7) 


The channel degree-of-freedom Va is related to the elastic 
enhancement factor Wa = 1 + 2/Va- 


3. KKM 

The model of Kawai, Kerman, McVoy [5] is very differ¬ 
ent from the MC approach of HRTW or Moldauer. The 
S'-matrix is expressed in terms of the optical S'-matrix 
background, in which the energy average of the reso¬ 
nance sum part will be zero. The optical model (or the 
coupled-channels optical model) yields Satchler’s trans¬ 
mission matrix [36], and a new hermitian matrix X in 
channel space is defined as 


includes the elastic enhancement and the width fluctua¬ 
tion correction at the same time, which is one of the rea¬ 
sons we defined the width fluctuation correction factor 
by Eq. namely the cross section ratio to the Hauser- 
Feshbach formula. 


5. Asymptotic expansion 

An asymptotic expansion of the GOE triple-integral 
formula in powers of 1/T is given by |34] 

{aab)^{l + 5ab)^A + 25ab%B, (A.IO) 

1 J ^ 

where 

A = 1 -I- — -I- —^ {T,2 — {Ta + Tb)} 

+ As2{S2-(r, + r,)} 

+ ^ (Tf + TaTb + Ti - Ea), (A. 11) 

S = (1 - T^a) |l - |(1 + 2Ta) + |S2|, (A.12) 

and 

^2 = ^ E S3 = ^ E Sc ■ (A.13) 

C C 

6. Ernebjerg and Herman 


X = -]^tvX + {{ivX/2f+ pY^‘^ . (A.8) 

In the overlapping resonance limit {T/D 3> 1), the av¬ 
erage cross section is written in terms of the A-matrix 
as 


Ernebjerg and Herman |22j generated a quasi-random 
set of transmission coefficients, and compared the sim¬ 
ulated cross sections with Eqs. (A.l), (A.4|, and (351. 
They obtained a new parameterization of the channel 
degree-of-freedom 


i^ab) — A^aa^bb 4 ” ^ab^ba 


(A.9) 


Since Eq. ( A.8| ) is a non-linear equation in A, one has 
to solve it by an iterative procedure (SSj. When {S) is 
diagonal (no direct channel), KKM yields an elastic en¬ 
hancement factor Wa = 2. In other words, KKM gives 
the correct asymptotic value in the Ericson regime. That 
same statement applies in the case of direct reactions. 
This is seen using the EW transformation. 


where 


1 

1 + f{Ta)T9(T.) ’ 


(A.I4) 


fiTa) 

9{Ta) 


0.177 

1 - 20.337ra ’ 
i + 3.i48r,(i-r,). 


(A.15) 

(A.16) 


7. Kawano and Talou 


4. GOE 


The analytical expression of the correct Hauser- 
Feshbach cross section, i.e. an analytical average over the 
GOE resonance parameter distributions, was given by 
Verbaarschot, Weidenmiiller, and Zirnbauer m. which 
is the so-called triple-integral of Eq. (351. The result 


Similar to Ernebjerg and Herman’s attempt, the GOE 
triple-integral calculation can be well-approximated 
by putting the following channel degree-of-freedom in 
Moldauer’s method 


Va = 2- 


1 

i+7 ’ 


/ = aPi 


Ta+T 


(A.17) 


1 -T, 
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They obtained 


obtained with 




/T-TaV 


a = 0.0288ra + 0.246 , 

(A.18) 

/ = 3a/?2Xa f ] , 

(A.20) 

/3i = l + 2.5r„(l-r„)exp(-2r) . 

(A.19) 

132 = 1 + 2.5Ta{l - To) exp(-4T) , 

(A.21) 

In the special case of T < 2Ta, a better fit 

can be 

where q = 0.4x°-^ and Xa = Ta/(1 — To). 




